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Abstract 



Coalescing neutron star binaries are believed to be the most reliable 
sources for ground-based detectors of gravitational waves and likely progen- 
itors of short gamma-ray bursts. In the process of coalescence, magnetic 
fields of neutron stars can induce interesting observational manifestations 
and affect the form of gravitational wave signal. In this papaer we use the 
population synthesis method to model the expected distribution of neutron 
star magnetic fields during the coalescence under different assumptions on 
the initial parameters of neutron stars and their magnetic field evolution. 
We discuss possible elecotrmagnetic phenomena preceding the coalescence 
of magnetized neutron star binaries and the effect of magnetic field on the 
gravitational wave signal. We find that a log-normal (Gaussian in loga- 
rithms) distribution of the initial magnetic fields of neutron stars, which 
agrees with observed properties of radio pulsars, produces the distribution 
of the magnetic field energy during the coalescence that adequately describes 
the observed luminosity function of short gamma-ray bursts under different 
assumptions on the field evolution and initial parameters of neutron stars. 
This agreement lends further support to the model of coalescing neutron star 
binaries as progenitors of gamma-ray bursts. 



1 Introduction 

Neutron stars (NS) are end products of thermonuclear evolution of stellar 
cores with a mass above ~ 8 — lOM©. The birth of a NS in the core collapse 
is accompanied by the ejection of the outer stellar envelope which is observed 
as type II or Ib/c supernova, depending on the evolutionalry status of the 
star immediately before the collapse. Young NS observed as radio pulsars 
demonstrate large space velocities (typically several hundreds km/s, see e.g. 
Hobbs et al. 2005), which is commonly explained by a hypothetical kick 
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velocity acquired by a nascent NS during the collapse. Different physical 
mechanisms of the kicks are duiscussed in the literature (e.g. Lai 2001). 

It is well known that astrophysical manifestations of a NS at all evolutionary 
stages are determined by its spin, magnetic field, and the interaction with 
the of the surrounding medium (see Lipunov 1992). The evolution of a NS in 
a binary system is even more complicated by the presence of the secondary 
component. So the properties of the galactic population of NS in binaries are 
subjected to uncertainties in the theory of evolution of binary stars, such as 
the initial orbital parameters of binary stars, uncertainties in the description 
of the evolutionary mass loss of normal stars via stellar winds, the efficiency 
of common envelops, etc. Binary NS observed as a specific subclass of binary 
pulsars are of special interest. Close binary NS (also called rclativistic double 
stars) are considered as best sources of gravitationla waves to ve observed by 
the present-day gravitational wave detectors (Grishchuk et al. 2001). The 
evolution of relativisitc binaries and binary pulsars is considereed in more 
detail in many reviews (see, e.g., Postnov and Yungelson (2006), Bisnovaty- 
Kogan (2006)). The coalescence of binary NS due to secular decrease of 
the orbit because of orbital angular momentum removal by gravitational 
waves can be the primary formation channel of short hard gamma-ray bursts 
(Bhnnikov et al. 1984, Eichler et al. 1989, see Nakar 2007 for extensive 
review) . 

Magnetic fields of coalescing NS binaries can lead to intersting observational 
manifestations preceding the gravitational wave burst (e.g. Lipunov and 
Panchenko 1996, Lipunova et al 1997), and can play a significant role in shap- 
ing the expected form of gravitational wave signal (Lipunova and Lipunov 
1998). The recent general rclativistic numerical calculations of coalescence 
of magnetized double NS (Yuk Tun Li et al. 2008) cofirmed the crucial role 
of strong magnetic field in delaying the collapse of an intermediate super- 
massive object that forms during the coalescence process. Obviously, the 
magnetic field cannot be neglected when calculating the gravitational-wave 
templates from colaescing NS binaries, which is very important for ongoing 
gravitational wave detectors (Braginsky 2008, Abbott et al. 2008). 

The population synthesis method is a powerful tool to study astrophysical 
manifestations of NSs. It enables calculation of large ensembles of objects and 
testing different scenarios of their evolution (see, e.g., Popov and Prokhorov 
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2007 for more detail). This method, with different degree of account for fea- 
tures of binary star evolution, has been successfully applied to study the evo- 
lution of coalescing NS binaries by differnt authors (Tutukov and Yungleson 
1994, Lipunov et al. 1997, Portegies Zwart and Yungelson 1998, Belczynski 
et al. 2008). In recent papers (Postnov and Kuranov 2008, Kuranov et al 
2009), the population synthesis method was apphed to study the effect of 
the possible kick-spin correlation on the coalescence rate of binary NSs and 
the observed correlation of the spin axis of radio pulsars with the direction 
of their space velocities. 

In the present work we use the Scenario Mashine code (see the detailed de- 
scription in Lipunov et al. 2009) to calculate the expected distribution of 
magnetic fields in coalescing binary NSs under different assumptions on the 
initial distribution of NS magnetic field and its evolution (no field decay or 
accretion-induced decay). The calculated distribution of magnetic fields is 
comapared with the observed luminosity function of short gamma-ray bursts. 
It is shown that for the initial log-normal distribution of NS magnetic fields, 
which well describes observational poroperties of radio pulsars, the calcu- 
lated distribution of the magnetic energy (both the energy of the poloidal 
magnetic field and the total magnetic field including the toroidal compo- 
nent) in coalescing binaries well fits the observed luminosity function of of 
short gamma-ray bursts if one assumes that the total energy release in a 
gamma-ray burst formed during the coalescence is proportional to the total 
magnetic field energy. The agreement is achieved in different models of the 
NS magnetic field evolution (the field decay time, accreiton-induced decay) 
and initial p[arameters of young NS (the amplitude of natal kick velocity and 
its possible correlation with the NS spin) . 



2 THE MODEL 
2. 1 Magnetic fields of neutron stars 

The initial magnetic field of neutron stars is one of the main parameters of 
our calculations. We used two possible forms of the initial magnetic field 
distribution of NSs: 1) log-normal distrubution of dipole magnetic moments 
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and 2) plain in logarithm distribution: 
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The first law directly follows from the analysis of the observed properties of 
radio pulsars and the field estimation using pulsar spin period and its first 
derivative (see the ATNF pulsar catalogue, Hobbs and Manchester 2009). 
The second law covers a broad range of magnetic field values and assumes 
an equal birth probability for strong-field NSs (like magnetars) and low- 
field NSs (like anti-magnet ars, i.e. radio pulsars with low surface fields and 
slow rotation). Indeed, some recent studies (e.g. Woods 2008) suggest the 
birthrate of magnetars comparable with that of ordinary pulsars, and one of 
the youngest NSs in the center of the supernova remnant Cas A is likely to be 
an anti-magnetar (Gotthelf and Halpern 2008). Of course, the NS magnetic 
filed at birth can depend on the NS parameters (for example, on its angular 
momentum) , but in view of the existing uncertainty in NS formation theory 
we do not consider such a complication of the initial conditions. 

In our calculations we adopted the following parameters of the initial dis- 
tribution of the NS dipole magnetic moment (Faucher-Giguere and Kaspi 



(the NS surface magnetic field B is calculated for the assumed NS radius 
Rm = 10 km). 



The evolution of the NS magnetic field is poorly known. It follows from 
general theoretical considerations that poloidal magnetic fields of single NSs 
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cr^ = 0.55, 
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2.2 Magnetic field evolution 
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must decay due to Ohmic losses in the crust (Flowers and Ruderman, 1977) 
with the characteristic time longer than ~ lO''' yrs. Observational tests of this 
hypothesis yield controversial results. For example, the analysis of statistical 
properties of radio pulsars leads some authors to conclude that the magnetic 
field of NSs must decay with the characteristic time ~ 3 mln yrs (Gontier et al 
2004), while another authors argue that statistical properties of radio pulsars 
do not require any magnetic field decay on timescales as long as comparable 
to the pulsar life time ~ 100 mln yrs (Faucher-Giguere and Kaspi 2006). 

In our calculations we used the following parametrization of the NS magnetic 
field evolution with time: 



The characterisitc decay time of the dipole magnetic field was chosen from 
T = 10'' yrs to oo (no decaY), and the minimum possible NS magnetic dipole 
moment was taken to be /Xmin = 10^^ G-cm^, which corresponds to a NS 
surface field of S ~ 10« G. 

Indeed, observations and theoretical arguments suggest that the surface poloidal 
magnetic field of an accreting NS in a close binary system must asymptoti- 
cally decay (Bisnovaty-Kogan and Romberg 1974, Taam and Van den Heuvel 
1986, Romani 1990, Bhattacharia 2002) down to a lower hmit of 10^ - 10^ G 
(Kulkarni 1986, Van den Heuvel et al. 1986). The magnetic field decay in 
accreting NSs was accounted for in the form (Romani 1990): 



where /jq and /j, are the value of the NS dipole magnetic moment before and 
after the accretion stage, respectively. Am is the mass of accreted matter (in 
solar masses). 



The stability analysis of an axisymmetric magnetic field (Prenderghast 1956, 
Tayler 1973, Wright 1973) suggests that a substantial fraction of the total 
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2.3 Toroidal field of NSs 
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magnetic field energy of a NS is contained in tfie toroidal component. Each 
field component (poloidal or toroidal) are separately unstable but can be 
stabilized by the presence of another component. The stability criterion of 
an axisymmetric magnetic field depends on several parameters and for a NS 
can be expressed in terms of the ratio of the energy of the poloidal component 
Ep to the total magnetic energy of the star = Ep + E^ (see Braithwaite 
2008): the upper limit for stability is {Ep/ E^) ~ 0.8; the lower limit for 
stability is {Ep/ Em)crit ~ '^O^Em/\Egr\. Here Egr is the gravitational binding 
energy (~ 10^^ erg). In our calculations we assumed that the NS magnetic 
field always satisfies this stability criterion, i.e. for a given value of Ep{t) the 
total magnetic field energy Ejn{t) falls within the stability limits 



The collapse asymmetry leads to additional kick velocity to a newborn NS. 
We modeled this velocity as follows: the direction of the kick velocity vector 
was taken to lie randomly within a cone with angle 9 around the rotational 
axis of the pre-supernova, and its absolute value was assumed to be dis- 
tributed according to a Maxwellian law: 



where Vm is the characteristic velocity. 

We have used two frequently considered kick velocity distributions (see the 
discussion in review Postnov and Yungelson (2006) and in the recent paper 
Martin et al. 2009). 

Model . Single-mode Maxwellian distribution. All NSs acquire kick velocity 
at birth. The parameter Vm varies within the range 100 — 500 km/s. 

Model . Bimodal Maxwellian distribution. For the components of close 
binary system (here we call a system as close binary if one of the components 



1.25Ep{t) < Em{t) < 10-^ ^Ep{t)\E, 



(6) 



2.4 Kick velocity during NS or BH formation 
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filled its Roche lobe and mass transfer occurred at some stages preceding the 
core collapse supernova of the primary component) with initial masses SM©- 
IIMq, the electron-capture supernova from degenerate 0-Ne-Mg stellar core 
is assumed. In that case the natal kick is small or zero. This assumptions 
allows one to explain some observational properties of X-ray Be-systems and 
binary pulsars (see Podsiadlowski et al 2004, Van den Heuvel 2004, 2007). 
In Model B with bimodal kick we assumed Vm = 30 / for a pre-supernova 
with the initial mass from the range 8Mq-11Mq and the mass transfer stage 
took place before the collapse. In other cases (i.e. for NS born from single 
progenitors or binary components with the initial mass above 11 solar masses) 
the parameter Vm varied within the range 100 — 500 km/s, as in Model A. 

500km/s, > IOM0, Model A 

500km/s, > IIM0, Model B 

30km/s, 8Mq-11Mq, closebinary 

For both kick models we assumed all single NSs and NSs in wide binaries to 
have the initial mass Mns = 1.4M0. In model B, NSs formed from e-capture 
SN are assumed to have smaller initial masses Mns = 1.25M0. 




3 RESULTS OF POPULATION SYNTHESIS CALCULATIONS 



The models of the natal kick and magnetic field evolution of netron stars are 
listed in Table 1. The results of calculations of different models are presented 
in Fig. 1-5. For each model, calculations were performed for isotropic kicks 
(corresponding to ^ = 90°) and kicks collimated with the rotational axis of 
the NS within the cone 9 = 10°. This value of 9 is suggested by the recent 
analysis of the spin- velocity alignment of radio pulsars (Kuranov et al. 2009). 

Fig. [T] shows the integral distribution of the binary NS coalescence rate as a 
function of the time of coalescence (determined as the time interval from the 
formation of a double NS binary until its merging due to gravitational wave 
emission) for different models of NS formation (models A and B). The total 
(over galactic lifetime) binary coalescence rate is a few times higher in Model 
B than in Model A (see also Postnov and Kuranov 2008), which is naturally 



7 



related to decrease in the binary disruption probability during formation of 
low-kick NSs in Model B. A strong dependence of the NS binary coalescence 
rate on the kick-spin alignment effect characterized by the value of the angle 9 
can also bee seen. For small 9 (strong alignment) the mean time of NS binary 
coalescence increases, so mainly old NSs coalesce. With increasing angle 
9 (the kick isotropisation) the fraction of young coalescing NSs with high 
magnetic field increases, since for isotropic kicks there is a certain probability 
for the NS to get kick directed oppositely to the orbital motion at the moment 
of explosion. Correspondingly, the fraction of coalescing highly magnetized 
NSs is maximum for isotropic kicks when 9 = 90°, see Fig. |2j 

Fig.[2]shows the distribution of NS dipole mometa (in units /iso = /i/(10'^° G cm^ 
at the coalescence for the tight kick-spin alignment (upper panel) and the 
isotropic kick (bottom panel) for a log-normal initial magnetic field distribu- 
tion (Model BG8). In Fig. [3] an integral distribution of the minimum (left 
plots) and maximum (right plots) magnetic moments of the components of 
coalescing double NSs are presented for different models. Plots in upper and 
bottom panels are calculated for the initial log-normal (Model BOB) and log- 
fiat (Model BF8) field distributions, respectively. The two-dimensional dis- 
tribution of magnetic moments of both components at coalescence is plotted 
in Fig. m for the log- normal initial distribution (model BF8) . The minimum 
E"™*" and maximum -E™"^ values of the total magnetic field energy of the co- 
alescing binary components (satisfied to the stability criterion (^) in Model 
BG8 is shown in Fig. [51 



4 DISCUSSION 



The calculated magnetic field distribution of coalescing NSs has several as- 
trophysical implications. 

1) Electromagnetic precursors of short gamma-ray bursts. 

A colaescing NS binary that leads to a powerful gamma-ray burst (Blinnikov 
et al. 1984, Eichler et al. 1989, see Nakar 2007 for a review) can be preceded 
by an electromagnetic burst (Lipunov and Panchenko 1996). Assuming the 
Goldreich- Julian current in the relativistic wind along open magnetic field 
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lines and the light cylinder radius Ri = cjVL ^ a, which are determined dy 
the rotation of the magnetized NS, we obtain the standard estimation of the 
total EM losses (e.g. Beskin 2009) 

Wem ~ iytR/cYB'^R^c ~ nVVc^ (8) 

(here Vl = 2it/P is the spin frequency of the NS, R is its radius, /i = BR^ 
is the NS dipole magnetic moment, a is the radius of the circular orbit, c is 
the speed of light). Assuming that the NS rotation is synchronized with its 
orbital relvolution at final stages before the coalescence, i.e. = 2GM/a^ 
(the 3-d Kepler law for two equal-mass NS in the circular orbit with semi- 
major axis a), we obtain 

VTem ~ 3 X lO39(erg/s)/i^o(M/M0)2(a/lOW)-6 . 

The characteristic values of the magnetic field of the second (younger) NS in 
the binary in the considered models is fi3o{NSII) ~ 10 — 10"'' (see Fig. [3]), 
and we obtain a total EM power emitted before the coalescence within the 
range from ~ 10^^ to ~ 10^^ erg/s, with the energy release strongly increasing 
with orbital decrease as 

In the relativistic particle wind flowing along open fleld lines plasma instabil- 
ity can develop to generate high-frequency radio emission, as in radio pulsars. 
The efficiency of conversion of the total power released into radio emission 
is fairly small in radio pulsars, ~ 10~^ — 10~^, so even during last orbital 
revolutions before the coalescence the expected power of a radio burst is 
hardly to be more than ~ 10^^ erg/s. In principle, this power is suffcient 
to explain properties of the unique extragalactic 5-ms radio burst (Lorimer 
et al 2007). However, such a powerful radio burst would require a significan 
dipole magnetic field of the NS at coalescence, fi^o ^ 1. Fig. Oshows that for 
the most favorable models with narrow kick-spin alignment the fraction of 
such highly magnetized NS is of order of a few per cents of the total number 
of coalescing NS binaries. This aggravates the low-event-rate problem for 
the explanation of the obsserved ms raio burst by the binary NS coalescence 
model (Lorimer et al 2007, Popov and Postnov 2007). On the other hand, 
in this model there is no problem for the highly coherent radio emission to 
go out of the generation region due to induced scattering (Lyubarsky 2008), 
which must be overcome in other models for this radio burst (e.g. Egorov 
and Postnov 2009). The fact that the observed ms radio burst was not asso- 
ciated with a gamma-ray burst could be due to different beaming angles for 
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the pre-outburst radio emission and the subsequent gamma-ray burst. Low 
statistics (one event), however, makes the statistical arguments very weak, so 
the possibihty to explain this radio burst by mechanism proposed by Lipunov 
and Panchenko (1996) still remains. 

2) The effect of the magnetic field of coalescing NSs on the form 
of the gravitational wave signal 

Yuk Tun Li et al (2008) studied numerically the effect of a strong magnetic 
field of the components of a coalescing NS binary on the form of the emitted 
gravitational wave (GW) signal at the stage preceding the merging of the 
components. In their numerical relativistic calculations, these authors used 
models of non-rotating NSs with a strong [B ~ 10^^ G) poloidal magnetic 
field. The difference in the GW signal form turned out to be significant 
if the final collapse of the supramassive differentially rotating NS (formed 
in the merging) into a black hole is magnetically delayed. Our calculations 
show that even for the initial high (magnetar-like) magnetic fields, the vast 
majority of coalescing NSs has significantly (by orders of magnitude) smaller 
magnetic fields if the field decay is assumed (Fig. H]). However, internal 
toroidal magnetic field can be much higher than the poloidal component, so 
that the energy of the magnetic field at coalescence can be mostly determined 
by the toroidal component and be as high as 10^^ ergs (and even more if 
one assumes the initial fields as high as ~ 10^^ — 10^^ G). Although such 
a field can not be too significant dynamically {Em/\Egr\ ~ 10~^ — 10~^), 
it may serve as a seed field for the magnetic field generation in a merged 
differentially rotating object. In principle, following arguments by Spruit 
(2008), the energy of differential rotation of the pre-collapsed merged object 
may entirely be transformed into the magnetic field energy, so that 

B^R^ ^ {An/nyErot, (9) 

where R is the radius of the object, the factor AQ/Q takes into account 
the differential rotation and Erot is the rotational energy. Assuming similar 
field generation processes to operate in a differentially rotating proto-NS and 
the pre-collapsing post-merging object, the magnetic field can be of order to 
Bmax,coal ~ Bmax,NS {,Erot,coai/Erot,NsY^'^- Limiting Erot,NS from above by the 
value ~ 10^^ ergs (Spruit 2008) and assuming the virial value of Erot,coai of 
order of the binding energy of the NS ~ 10^'^ ergs, we obtain the potentially 
possible energy of the poloidal magnetic field of pre-collapsing post-merging 
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object about ~ 10 times higher than for a young NS. This energy can be 
dynamically significant for the subsequent evolution of this object and must 
be taken into account in calculations of the expected form of the GW signal 
at post-merging phases. 

3) The luminosity function of short gamma-ray bursts 

Probably, the most intriguing implication of our calculations relates to short 
gamma-ray bursts. Presently, the merging of compact NS binaries is thought 
to be the principal formation channel of short hard gamma-ray bursts (Blin- 
nikov et al. 1984, Eichler et al. 1989, see Nakar 2007 for a recent review). The 
plausible physical mechanism of generation of hard electromagnetic gamma- 
radiation during the NS binary coalescence is related to the formation of 
ultra-relativisitc MHD jet (Usov 1992, see recent numerical calculations by 
Komissarov et al. 2009). The NS magnetic field can be an important initial 
parameter for such a jet formation. Indeed, the power of energy release in the 
MHD-mechanism is roughly L ~ B'^uj^R^, where u) is the characteristic spin 
frequency, R is the radius and B is the magnetic field strength. For different 
coalescing NS binaries parameters u and R are most likely similar, while the 
value of the magnetic field strength can be significantly different. Of course, 
as we argued above, the magnetic field can strongly increase during the co- 
alescence process, however immediately before the merging NS fields serve 
as the initial conditions. Thus, there are physical grounds to expect that the 
distribution of luminositites of short gamma-ray bursts ^{L) will reflect the 
distribution of magnetic fields of coalescing NSs: d^{L)/dL oc dN/dB. 

We shall consider two different possibilities. 

1) $(L) reflects distribution of the energy of the NS poloidal fleld Ep Bp 
before the coalescence, 

$(L) oc dN/dEp oc dN/d{Bl) . 

and 

2) $(-L) reflects the distribution of the total magnetic fleld energy (including 
the toroidal component) Em = Ep + Et. According to the stability criterion 
for the poloidal magnetic fleld ([6]), the maximum total energy of the magnetic 
fleld scales as E^ ~ \/ \ Egr \ Ep oc Bp, i.e. for a given value of the poloidal fleld 
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strength Bp, the distribution of the maximum possible magnetic field energy 
in a stable configuration is proportional to that of the poloidal magnetic field 
strength, 

$(L) oc dN/dEm oc dN/dBp . 

Our calculations of magnetic fields of coalescing NS binaries allows testing 
these hypotheses. To this aim, we shall use the luminosity function of short 
gamma-ray bursts calculated using observational data for short gamma-ray 
bursts obtained by space observatories Swift, INTEGRAL, Fermi (Kann et 
al 2008). In Figs. [6] -[7] the integral luminosity funciton of short gamma-ray 
bursts A^(> L) is shown by the solid black line connecting solid diamonds. 
The dashed lines show 1-a errors of the total GRB luminosity estimation 
(assuming spherical symmetry of emission) taken from Table 1 of Kann et al. 
(2008). In Fig. [6] we compare the hypothesis oc for different initial NS 
parameters (the natal kick velcoity parameter V^, spin- kick alignment angle 
9) in our Models AG and BG (in which the NS magnetic field can decay 
only during accretion) with log-normal initial distribution of the poloidal 
magnetic field. The upper and lower panel of Fig. show the comparison 
with the poloidal component of the second (younger) NS in the binary, the 
middle panel show the comparison with the field of the first (older) NS. It 
is seen that the observed short GRB luminosity distribution is in agreement 
with the distribution of magnetic fields of NS before the coalescence for any 
assumptions on the value and the direction of the natal NS kick velocity. 

In Fig. [7] we compare the observed luminosity function of short gamma-ray 
bursts with the distribution of the total energy of the magnetic field before the 
coalescence (i.e. the hypothesis oc E^ oc Bp) for models with exponential 
field decay AG8 (the upper panel) and BG8 (the bottom panel), as well as 
for a log-fiat initial field distribution. For these models, the agreement with 
observations is reached for all values of narrow-collimated kicks [9 = 10^), 
and cannot be reached at high values of the isotropic kicks (Vm = 500, 300 
km/s and 9 = 90°). The parameter of the Kolmogorov-Smirnov test for 
different models is listed in Table 2. 

Note that in models with a logarithmically flat initial magnetic fleld distri- 
bution (our models AF, AF8, BF, BF8) the observed luminosity function of 
short gamma-ray bursts can not be described for any value and direction of 
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the NS kick. 



5 CONCLUSIONS 



Using the population synthesis method (the Scenario Mashine code), we have 
calculated the expected distribution of magnetic fields of coalescing binary 
neutron stars for various initial distributions and models of NS formation. 
We have taken into account both poloidal and toroidal components of the 
field satisfying the field stabihty criterion. The most important parame- 
ters of NS formation include the form (single-mode, bimodal Maxellian dis- 
tribuitons) and amplitude (30-500 km/s) of the natal kick velocity and its 
possible alignment with the NS spin. We have taken into account the NS 
magnetic field evolution (no decay, exponential decay until a bottom value 
of 10^ G, accretion-induced decay). We discuss possible observational man- 
ifestations of the strong magnetic field of coalescing neutron stars (the elec- 
tromagnetic precursors) and the effect of the magnetic field on the form of 
gravitational wave signal. 

We have shown that the observed luminosity function of short/hard gamma- 
ray bursts is very well fitted by the magnetic field distribution of the coalesc- 
ing neutron star binaries (the Kolmogorov-Smirnov test of order of one) for 
different initial NS formation models if the initial magneitc field distribution 
is taken in the log-normal form with parameters inferred from the radio pulsar 
observations. The obtained agreement of the distribution of magnetic fields 
of coalescing NS with the observed luminosity funciton of short gamma-ray 
bursts lends further support to the pioneer idea by Blinnikov et al. (1984) 
that coalescing NS binaries can be progenitors of cosmic gamma-ray bursts. 



The work is supprted by RFBR grant 07-02-00961. 
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Table 1. Models of the intital magnetic fields and natal kicks of NS 





Initial distribution 


Bp field decay 


Via Initial mass 




Model 


of magnetic field 


time, yrs 


km s"-*^ 


of NS progenitor 
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log-normal 
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Table 2. Kolmogorov-Smirnov test for the luminosity function of short gamma-ray 

bursts 



Tested 



Model 


liy])()triesis 





i;n=100 km s- 


^ \\,,=300 km s- 


^ \\,,=rm km 


AG 


Lry (X Ep 


10° 


0.904 


0.995 


0.525 
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BG 


Lry OC Ep 
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0.997 
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Lry OC Ep 
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0.893 


0.947 
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Lj oc E rn 


10° 
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0.999 
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0. 169 


2.6e-3 


3.8e-l 


BG8 


L-y oc Em 


10° 


1.000 


1.000 


1.000 




L-y oc Era 


90° 


0.999 


0.319 


0.071 
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Figure 1. Integral distribution function of the binary NS coalescence rate as a func- 
tion of time between the birth and merging for different NS formation parameters. 
The rate is normalized to the galactic star formation rate 3Mq Upper plots: 
model A, bottom plots: model B. Plots to the left: the natal kick-spin alignment 
within the cone with angle 9 = 10°. Plots to the right: the isotropical natal kick 
{9 = 90°). 
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Figure 2. The distribution function of dipole magnetic moments of coalescing 
meutron stars (NSI denotes the first (older) NS, NS II - the second (younger) NS 
in the binary) by the moment of coalescence in Model BG8. The upper plot is 
calcTilatcd for narrow collimated kicks within 9 = 10°. The bottom plots shows 
the result for isotropic kicks {9 = 90°). 
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Figure 3. The integral distribution function of the minimum (to the left) and 
maximum (to the right) dipole moment of two binary NS components by the time 
of coalescence. The upper and bottom plots are calculated for Model BG8 and 
BF8, respectively. The angle 9 for natal kick-spin correlation is also shown. 




Figure 4. The distribution function of the dipole magnetic moments of the poloidal 
magnetic field and the coalescence time in Model BG8. The upper and bottom 
plot shows the distribution for the first (older) NSI and the second (younger) NSII, 
respectively. Left and right plots are for narrow-collimated {6 = 10°) and isotropic 
{9 = 90°) natal kicks, respectively. 
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Figure 5. The minimum (-B^*", the upper panel) and maximum (E^^ ^ the bottom 
panel) total magnetic field energy of coalescing binary NS components by the time 
of coalescence. Model BG8. 



23 




49 50 51 52 53 54 




lg(L/erg) 



Figure 6. The comparison of the hypothesis oc for the luminosity of short 
gamma-ray bursts for different NS formation parameters in models AG and BG. 
On the top and bottom panel the comparison with the distribution of the poloidal 
field of the second (younger) NS in the binary is shown, while the plot on the middle 
panel is calculated for the field distribution of the first (older) NS. The observed 
integral distribution function of short gamma-ray bursts A^(> L) as computed 
from data presented in Table 1 of Kann et al. (2008) is shown by the solid line 
connecting the filled diamonds. The dashed lines indicate l-a errors of the gamma- 
rya burst luminosity L-y. 
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Figure 7. The comparison of the hypothesis oc E„i for the luminosity of short 
gamma-ray bursts in models with exponential field decay AGS (the upper panel) 
and BG8 (the bottom panel) for the initial log-normal distribution of the poloidal 
magnetic field component of NS. 
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